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Thermoelectric  materials  interconvert  thermal  gradients  and 
electric  fields  for  power  generation  or  for  refrigeration1,2. 
Thermoelectrics  currently  find  only  niche  applications  because 
of  their  limited  efficiency,  which  is  measured  by  the  dimensionless 
parameter  ZT — a  function  of  the  Seebeck  coefficient  or  ther¬ 
moelectric  power,  and  of  the  electrical  and  thermal  conductivities. 
Maximizing  ZT  is  challenging  because  optimizing  one  physical 
parameter  often  adversely  affects  another3.  Several  groups  have 
achieved  significant  improvements  in  ZT  through  multi-component 
nanostructured  thermoelectrics4-6,  such  as  Bi2Te3/Sb2Te3  thin-film 
superlattices,  or  embedded  PbSeTe  quantum  dot  superlattices. 
Here  we  report  efficient  thermoelectric  performance  from  the 
single-component  system  of  silicon  nanowires  for  cross-sectional 
areas  of  10  nm  X  20  nm  and  20  nm  X  20  nm.  By  varying  the  nano¬ 
wire  size  and  impurity  doping  levels,  ZT  values  representing  an 
approximately  100-fold  improvement  over  bulk  Si  are  achieved 
over  a  broad  temperature  range,  including  ZT  ~  1  at  200  K. 
Independent  measurements  of  the  Seebeck  coefficient,  the  elec¬ 
trical  conductivity  and  the  thermal  conductivity,  combined  with 
theory,  indicate  that  the  improved  efficiency  originates  from 
phonon  effects.  These  results  are  expected  to  apply  to  other  classes 
of  semiconductor  nanomaterials. 

The  most  efficient  thermoelectrics  have  historically  been  heavily 
doped  semiconductors  because  the  Pauli  principle  restricts  the  heat¬ 
carrying  electrons  to  be  close  to  the  Fermi  energy1  for  metals.  The 
Wiedemann-Franz  law,  kJoT=  tc2/3 (fc/e)2  =  (156  |iVK-1)2,  where 
Ke  is  the  electronic  contribution  to  k,  constrains  ZT  =  S2cfTIk,  where 
S  is  the  Seebeck  coefficient  (or  thermoelectric  power,  measured  in 
VK  1),  and  a  and  K  are  the  electrical  and  thermal  conductivities, 
respectively.  Semiconductors  have  a  lower  density  of  carriers,  leading 
to  larger  S  values  and  a  k  value  that  is  dominated  by  phonons  (fcph), 
implying  that  the  electrical  and  thermal  conductivities  are  somewhat 
decoupled1,  k  can  be  reduced  by  using  bulk  semiconductors  of  high 
atomic  weight,  which  decreases  the  speed  of  sound.  However,  this 
strategy  has  not  yet  produced  materials  with  ZT>  1. 

For  a  metal  or  highly  doped  semiconductor,  S  is  proportional  to 
the  energy  derivative  of  the  density  of  electronic  states.  In  low-dimen¬ 
sional  (nanostructured)  systems  the  density  of  electronic  states  has 
sharp  peaks7-9  and,  theoretically,  a  high  thermopower.  Harnessing 
this  electronic  effect  to  produce  high-ZT  materials  has  had  only 
limited  success10,11.  However,  optimization  of  the  phonon  dynamics 
and  heat  transport  physics  in  nanostructured  systems  has  yielded 
results4-6.  Nanostructures  maybe  prepared  with  one  or  more  dimen¬ 
sions  smaller  than  the  mean  free  path  of  the  phonons  and  yet  larger 
than  that  of  electrons  and  holes.  This  potentially  reduces  k  without 
decreasing  a  (ref.  12).  Bulk  silicon  (Si)  is  a  poor  thermoelectric 
(■ZT300K  ~  0  .01;  ref.  13),  and  this  phonon  physics  is  important  for 
our  Si  nanowires,  in  which  the  electronic  structure  remains  bulk-like. 


Figure  1  shows  images  of  the  devices  and  the  Si  nanowires  we  used 
for  these  experiments.  Details  related  to  the  fabrication,  calibration 
and  experimental  measurements  are  presented  in  the  Supplementary 


Figure  1  |  Scanning  electron  micrographs  of  the  device  used  to  quantitate 
the  thermopower  and  electrical  and  thermal  conductivity  of  Si  nanowire 
arrays,  a,  This  false-colour  image  of  a  suspended  platform  shows  all 
electrical  connections.  The  central  green  area  is  the  Si  nanowire  array,  which 
is  not  resolved  at  this  magnification.  The  four-lead  yellow  electrodes  are  used 
for  thermometry  to  quantify  the  temperature  difference  across  the  nanowire 
array.  The  thermal  gradient  is  established  with  either  of  the  two  Joule  heaters 
(the  right-hand  heater  is  coloured  red).  The  yellow  and  blue  electrodes  are 
combined  to  carry  out  four-point  electrical  conductivity  measurements  on 
the  nanowires.  The  grey  region  underlying  the  nanowires  and  the  electrodes 
is  the  150-nm-thick  Si02  insulator  that  is  sandwiched  between  the  top 
Si(100)  single-crystal  film  from  which  the  nanowires  are  fabricated,  and  the 
underlying  Si  wafer.  The  underlying  Si  wafer  has  been  etched  back  to 
suspend  the  measurement  platform,  placing  the  background  of  this  image 
out  of  focus,  b,  Low-resolution  micrograph  of  the  suspended  platform.  The 
electrical  connections  radiate  outwards  and  support  the  device,  c,  High- 
resolution  image  of  an  array  of  20-nm-wide  Si  nanowires  with  a  Pt  electrode. 
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Figure  2  |  Factors  contributing  to  ZT  for  various  Si  nanowires.  All 
nanowires  are  20  nm  in  height,  a,  The  temperature  dependence  of  the 
thermal  conductivity  K,  presented  as  /Cbuli/t'nanowirc  to  highlight  the 
improvement  that  the  reduction  of  K  in  nanowires  lends  to  ZT.  Jtbuik  values, 
which  are  slightly  below  the  true  bulk  value  for  Si,  are  taken  from  an 
identically  measured  520  nm  X  35  nm-sized  film.  The  inset  scanning 
electron  microscope  micrographs  show  the  region  of  the  device  containing 
the  nanowires  before  (top)  and  after  (bottom)  the  XeF2  etch  to  remove  the 


nanowires,  b,  The  temperature  dependence  of  S2  for  20-nm-wide  Si 
nanowires  at  various  p-type  doping  concentrations  (indicated  on  the  graph) . 
Note  that  the  most  highly  doped  nanowires  (pink  line)  yield  a  thermopower 
similar  to  that  of  bulk  Si  doped  at  a  lower  level.  For  nanowires  doped  at 
slightly  higher  and  slightly  lower  concentrations  than  the  bulk,  S  peaks  near 
200  K.  This  is  a  consequence  of  the  one-dimensional  nature  of  the  Si 
nanowires. 


Information.  The  platform  permits  four-point  measurements  of  the 
electrical  conductivity  of  the  nanowires,  Joule  heating  to  establish  a 
thermal  gradient  across  the  nanowires,  and  four-point  thermometry 
to  quantify  that  gradient10,14.  The  resistance  of  the  four-point 
thermometry  electrodes  is  typically  two  orders  of  magnitude  smaller 
than  the  resistance  of  the  nanowires.  The  measurement  platform  is 
suspended  in  vacuum  to  allow  measurement  of  nanowire  thermal 
conductivity15,16.  For  all  measurements,  the  Si  nanowires  could  be 
selectively  removed  using  a  XeF2  etch,  thus  allowing  for  measure¬ 
ments  of  the  contributions  to  c,  S  and  K  from  the  platform  and  oxide 
substrate. 

There  are  several  ways  to  prepare  Si  nanowires,  including  materials 
methods  for  bulk  production17.  We  wanted  Si  nanowires  in  which  the 
dimensions,  impurity  doping  levels,  crystallographic  nature,  and  so 
on,  were  all  quantifiable  and  precisely  controlled.  We  used  the  super¬ 
lattice  nanowire  pattern  transfer  (SNAP)  method18,  which  translates 
the  atomic  control  over  the  layer  thickness  of  a  superlattice  into 
control  over  the  width  and  spacing  of  nano  wires.  Si  nanowires  made 
via  SNAP  inherit  their  impurity  dopant  concentrations  directly  from 
the  single-crystal  Si  epilayers  of  the  silicon-on-insulator  substrates 
from  which  they  are  fabricated19.  These  epilayers  were  20-  or  35-nm- 
thick  Si(100)  films  on  150  nm  of  Si02,  and  were  p-type  impurity 
(boron)  doped  using  diffusion-based  doping19.  Four-point  probe 
conductivity  measurements  of  the  silicon-on-insulator  films  were 
used  to  extract  dopant  concentrations.  We  prepared  nanowire  arrays 
several  micrometres  long,  with  lateral  width  X  thickness  dimensions 
of  10nmX20nm,  20nmX20nm  and  520nmX35nm.  The  last 
approximates  the  bulk  and,  in  fact,  measurements  on  the  sample 
obtained  bulk  values  for  S,  a  and  k.  Measurements  of  k  for  our 
nanowires  were  consistent  with  literature  values  for  materials  grown 
(round)  Si  nanowires16. 

All  values  of  S,  a  and  K  reported  here  are  normalized  to  individual 
nanowires,  although  each  experiment  used  a  known  number  of 
nanowires  ranging  from  10  to  400.  The  Si  microwires  and  nanowires 
were  prepared  using  the  same  substrates,  doping  methods,  and  so  on, 
but  different  patterning  methods  (electron-beam  lithography  versus 
SNAP). 

Measurements  of  k  and  S2  for  Si  nanowires  (and  microwires)  for 
different  nano  wire  sizes  and  doping  levels  are  presented  in  Fig.  2. 
More  complete  data  sets,  electrical  conductivity  data  and  a  statistical 
analysis  are  presented  in  the  Supplementary  Information.  The  nano¬ 
wire  electrical  conductivity  is  between  10  and  90%  of  the  bulk, 
depending  on  nanowire  dimensions.  A  reduced  a  probably  arises 
from  surface  scattering  of  charge  carriers19.  Nevertheless,  all  nano¬ 
wires  are  highly  doped  and  most  exhibit  metallic-like  conductivity 
(increasing  a  with  decreasing  T). 


The  temperature  dependence  of  k  for  a  micro  wire  and  10-  and 
20-nm-wide  nanowires  were  recorded  at  modest  statistical  resolution 
to  establish  trends.  This  data  indicated  that  K  drops  sharply  with 
shrinking  nanowire  cross-section  (Fig.  2a)  and  that  the  10-nm-wide 
nanowires  exhibited  a  k  value  (0.76±0.15Wm  *K  ')  that  was 
below  the  theoretical  limit  of  0.99  W  m  'K  1  for  bulk  Si  (ref.  20). 
Thus,  very  large  data  sets  were  collected  for  10-  and  20-nm-wide 
nanowires  to  allow  for  a  more  precise  determination  of  k. 

Our  observed  high  ZT  for  Si  nanowires  (Fig.  3)  occurs  because  K  is 
sharply  reduced  and  the  phonon  drag  component  of  the  thermo¬ 
power  Sph  becomes  large.  Below,  we  show  that  Sph  increases  because 
of  a  three-dimensional  to  one-dimensional  crossover  of  the  phonons 
participating  in  phonon  drag  and  decreasing  k.  However,  we  first 
discuss  why  our  measured  K  at  300  K  for  10-nm-wide  Si  nano  wires  is 
less  than  icm;n  (ref.  20). 

The  derivation  of  Kmin  assumes  that  the  minimum  path  length  of 
wavelength  2  phonons  is  2/2  and  that  the  phonons  are  described  by 
the  Debye  model  using  bulk  sound  speeds  with  no  optical  modes.  The 
2/2  value  is  an  order-of-magnitude  estimate  and  is  difficult  to  deter¬ 
mine  precisely,  much  like  the  minimum  electron  mean  free  path  used 
to  calculate  the  Mott-Ioffe-Regel  erm;rl.  Also,  Kmin  is  proportional  to 
the  transverse  and  longitudinal  acoustic  speeds  of  sound20.  These  are 
reduced  in  our  nanowires  at  long  wavelengths  because  the  modes 


Figure  3  |  Temperature  dependence  of  ZT  for  two  different  groups  of 
nanowires.  The  cross-sectional  area  of  the  nanowires,  and  the  p-type  doping 
level,  are  given.  The  20-nm-wide  nanowires  have  a  thermopower  that  is 
dominated  by  phonon  contributions,  and  a  ZT  value  —  1  is  achieved  near 
200  K.  The  smaller  (10-nm-wide)  nanowires  have  a  thermopower  that  is 
dominated  by  electronic  contributions.  The  ZT  at  350  K  is  calculated  using 
the  thermal  conductivity  value  for  the  10-nm-wide  nanowires  at  300  K.  The 
error  bars  represent  95%  confidence  limits. 
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become  one- dimensional.  The  ratio  of  the  one-dimensional  to  two- 
dimensional  longitudinal  speeds  of  sound  is  [(1  +  v)(l  —  2v)/ 
(1  —  v)]1/2  =  0.87  where  v  =  0.29  is  the  Poisson  ratio  of  Si.  The  trans¬ 
verse  acoustic  speed  goes  to  zero  at  long  wavelength  because  to  oc  k2  d 
where  d  is  the  nano  wire  width21.  Therefore,  the  bulk  icmin  estimate 
above  is  invalid  for  our  nanowires  and  values  smaller  than  Km;n  are 
attainable. 

For  all  but  the  most  highly  doped  nanowires,  S  peaks  near  200  K 
(Fig.  2b).  This  peak  is  unexpected:  similarly  doped  bulk  Si  exhibits  a 
gradual  decrease  in  S  as  T is  reduced  (green  trace).  For  T <  100  K,  a 
peaked  S(  T)  is  observed  for  metals  and  lightly  doped  semiconductors 
and  is  due  to  phonon  drag13,22,23. 

Phonon  drag  is  generally  assumed  to  vanish  with  decreasing  sam¬ 
ple  dimensions  because  the  phonon  path  length  is  limited  by  the 
sample  size24-26.  This  seems  to  eliminate  phonon  drag  as  the  reason 
for  the  peak  in  our  nano  wires.  We  show  below  that  the  phonon 
wavelengths  participating  in  drag  are  of  the  order  of  or  larger  than 
the  wire  width.  This  leads  to  a  three-dimensional  to  one-dimensional 
crossover  of  these  modes  and  removes  the  cross-sectional  wire 
dimensions  from  limiting  the  phonon  mean  path  (see  Fig.  4  inset). 
The  nanowire  boundaries  are  incorporated  into  the  one-dimensional 
mode  and  are  not  an  obstacle  to  phonon  propagation.  Therefore,  the 
limiting  size  becomes  the  wire  length  (~1  pm)  and  phonon  drag 
‘reappears’  at  very  small  dimensions. 

In  addition,  classical  elasticity  theory21  is  valid  for  the  phonon 
wavelengths  considered  here27,  leading  to  thermoelastic  damp¬ 
ing21,28'29  of  sound  waves  proportional  to  k.  Thus  Sph  is  further 
enhanced,  owing  to  the  observed  reduced  thermal  conductivity  k. 
A  detailed  discussion  is  in  the  Supplementary  Information. 

It  might  seem  that  elasticity  theory  leads  to  a  contradiction  because 
k  is  proportional  to  the  mean  phonon  lifetime.  If  the  phonon  life¬ 
times  increase  as  stated  above,  then  k  should  also  increase.  But 


Figure  4  |  Thermopower  calculation  plotted  along  with  experimental  data 
(black  points)  from  a  20-nm-wide  Si  nanowire  p-type  doped  at 
3  x  1019cm-3.  The  black  curve  is  the  fitted  expression  for  the  total 
thermopower  Se  +  Sph.  The  red  curve  is  the  phonon  contribution  Sph  and  the 
blue  line  is  the  electronic  term  Se  arising  from  the  fit.  The  fit  has  maximum 
error  6.1  pV  K  1  and  root-mean-square  error  1.8  pV  K  .  The  experimental 
error  bars  represent  95%  confidence  limits  and  at  150,  200  and  225  K  are 
smaller  than  the  data  points.  The  blue  data  points  are  experimental  values 
for  bulk  wires  (doping  2  X  1020  cm-3;  crosses),  10-nm-wide  nanowires 
(doping  7  X  1019  cm  3;  diamonds),  and  20-nm-wide  wires  (doping 
1.3  X  1020  cm-3;  triangles)  where  only  a  linear-T electronic  contribution  was 
found.  This  data  are  close  to  the  extracted  electronic  contribution  from  the 
black  data  points  (blue  line)  and  shows  that  the  fitted  linear  term  is 
reasonable.  The  drop  in  S  to  0  as  T— » 0  occurs  because  the  phonon  mean  free 
path  reaches  the  sample  size  and  the  specific  heat  tends  to  0  according  to  the 
third  law  of  thermodynamics.  The  inset  shows  the  character  of  a  three- 
dimensional  bulk  longitudinal  acoustic  phonon  mode  (top)  and  a  one¬ 
dimensional  mode  when  the  wavelength  is  larger  or  of  the  order  of  the  width 
(bottom).  The  one-dimensional  mode  incorporates  the  existence  of  the 
boundary  by  transverse  expansion  (or  compression)  for  longitudinal 
compression  (or  expansion).  The  ratio  of  the  transverse  strain  to  the 
longitudinal  strain  is  the  Poisson  ratio  (0.29  for  Si). 


because  the  elasticity  expression  is  only  valid  for  long-wavelength 
modes,  and  k  is  the  average  of  all  modes,  there  is  no  contradiction. 

We  now  consider  separately  the  electronic  and  phonon  contribu¬ 
tions  to  the  thermopower — S=  Se+Sph — for  the  nanowire  data  at 
T>200K.  Charge  carriers  dissipate  heat  to  the  lattice  through  a 
process  that  first  involves  momentum  conserving  (non-dissipative) 
electron-phonon  collisions.  The  phonons  that  contribute  to  phonon 
drag  cannot  have  a  wavelength  shorter  than  2min,  which  is  deter¬ 
mined  by  the  size  of  the  Fermi  surface.  Phonon  drag  is  observed  in 
metals  only  at  low  T  because  the  Fermi  surface  is  large  and  the  heat 
carrying  short-wavelength  phonons  have  short  lifetimes.  At  low  T 
(<20K),  Sph  oc  T3  from  the  phonon  specific  heat  (oc  r3).  For 
kT3>  ©Dj-bye,  the  specific  heat  becomes  constant  and  the  number 
of  phonons  available  for  phonon-phonon  scattering  is  oc  T,  leading 
to  Sph  oc  1/T  (ref.  1). 

For  p-type  Si,  the  holes  are  near  the  valence  band  maximum.  The 
phonon  drag  modes  are  acoustic  with  the  largest  wavevector 
kph  =  2kFermi  =  0.2  A  1  (for  impurity  doping  of  3X1019cm  3). 
The  shortest  wavelength  is  2ph  =  2nlk  =  31  A. 

Umklapp  (non-momentum-conserving)  phonon-phonon  scat¬ 
tering  processes  determine  the  rate  of  phonon  heat  dissipation. 
The  Debye  energy  0D  sets  the  energy  scale  for  Umklapp  scattering. 
The  number  of  Umklapp  phonons  available  to  dissipate  the  long- 
wavelength  phonons  is  given  by  the  Bose-Einstein  function: 


Nu=- 


1 


e0o/T  _  l 

leading  to  a  scattering  rate  l/rph  oc  Nv.  When  T3>  0D,  l/rph  oc  T. 
Because  0D  =  640  K  for  Si,  the  full  Bose-Einstein  expression  must  be 
applied  for  T £  350  K.  The  electronic  contribution  Se  is  estimated 
from  the  Mott  formula1: 

n2  k2T  S' din  g(£)\  . 

Se(T)=  M  =(283 pV  K-'XfcT/Ep) 

where  the  conductivity  derivative  equals  the  reciprocal  of  the  energy 
scale  over  which  it  varies  (the  Fermi  energy  EF).  Assuming  hole 
doping  occurs  in  the  heavier  Si  valence  band  (mass  0.49),  this  leads 
to  EF  =  0.072  eV  =  833  K  and  kF  =  0.1  A  1  for  the  number  density  of 
boron  dopant  atoms  n  =  3  X  1019cm  3.  Thus  S^(T)  =  aT  where 
a  =  0.34  pVK~2. 

The  T>  200  K  thermopower  data  of  the  20-nm-wide  wire  (doping 
n  =  3  X  1019  cm  3)  fits: 

S=  Se  +  Sph  =  aT+  fe[exp(@D/T)  —  1] 

where  a,  b  and  0D  are  varied  to  obtain  the  best  fit  (Fig.  4).  The  coeffi¬ 
cients  are  a  =  0.337  pVKT2,  b  -  22.1  pVK-1,  and  0D  =  534  K.  The 
coefficient  a  is  almost  identical  to  our  estimate  of  0.34  pV  K  2.  Thus 
phonon  drag  explains  the  observed  thermopower.  Consistent  with 
measurements24  of  phonon  drag  in  bulk  Si,  S  in  our  nanowires 
increases  significantly  at  lighter  doping.  This  data,  plus  a  fit  of  the 
Fig.  4  data  using  the  experimental  0D  =  640  K,  is  presented  in  the 
Supplementary  Information.  Rather  than  fitting  by  varying  a,  b  and 
0D,  we  fixed  0D  to  its  known  experimental  value  (640  K)  and  allowed 
only  a  and  b  to  vary. 

The  phonon  drag  contribution  to  S  is  of  the  form23,30: 

Sph®1 1 


V-T. 


Tph,  the  phonon  lifetime,  is  oc  1  Ik  from  elasticity  theory.  /(  is  the 
electron  mobility.  ZT scales  as  (neglecting  Se): 


^ph 


[.iTk 


ffocn/i 


ZT  oc  ■ 


/(Tk3 
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leading  to  increased  ZT with  decreasing  mobility.  This  is  the  opposite 
conclusion  reached  from  when  we  consider  only  Se  (ref.  7). 

We  have  combined  experiment  and  theory  to  demonstrate  that 
semiconductor  nanowires  can  be  designed  to  achieve  extremely  large 
enhancements  in  thermoelectric  efficiency,  and  we  have  shown  that 
the  temperature  of  maximum  efficiency  may  be  tuned  by  changing 
the  doping  and  the  nanowire  size.  Theory  indicates  that  similar 
improvements  should  be  achievable  for  other  semiconductor  nano¬ 
wire  systems  because  of  phonon  effects.  These  nano  wire  thermoelec¬ 
trics  may  find  applications  related  to  on-chip  heat  recovery,  cooling 
and  power  generation.  Additional  improvements  through  further 
optimization  of  nanowire  size,  doping  and  composition  should  be 
possible. 

METHODS  SUMMARY 

Single- crystalline  Si  nanowires  were  fabricated  using  the  SNAP  process18.  The 
nanowires  were  doped  p-type  using  a  boron-containing  spin-on  dopant  (see 
Supplementary  Information  for  details)19.  Electron-beam  lithography  was  used 
to  create  Ti/Pt  electrodes  for  the  electrical  and  heat  transport  measurements.  The 
entire  device  was  suspended  using  a  XeF2  etch,  leaving  the  nanowires  anchored  to 
a  thin  SiO  island  (see  Supplementary  Information  for  details).  The  nanowire 
electrical  conductivity  was  measured  by  a  Keithley  2400  using  a  four-point 
measurement  to  eliminate  contact  resistance.  For  measurement  of  S  and  k,  a 
temperature  difference  was  created  across  the  ends  of  the  nanowires  by  sourcing 
a  direct  current  through  one  of  the  resistive  heaters.  The  resistance  rise  of  each 
thermometer  was  recorded  simultaneously  using  a  lock-in  measurement 
(Stanford  Research  Systems  SRS-830)  as  the  temperature  was  ramped  upwards. 
The  resistance  of  the  thermometers  was  typically  two  orders  of  magnitude  smal¬ 
ler  than  the  nanowire  array.  For  measurement  of  S,  the  thermoelectric  voltage,  as 
a  response  to  the  temperature  difference,  was  recorded  using  a  Keithley  2182A 
nanovoltmeter.  A  difference  measurement  was  used  to  determine  k,  whereby  the 
k  value  of  the  nano  wires  plus  the  oxide  island  was  subtracted  from  the  k  value  of 
the  oxide  island.  The  thermal  conductivity  of  the  oxide  island  was  determined  by 
removing  the  nanowires  with  a  highly  selective  XeF2  etch  (see  Supplementary 
Information  for  details). 
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